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Many approximations witlriir deirsity-functional theory spuriously predict that a many-electron 
system can dissociate into fractionally charged fragments. Here, we revisit the case of dissociated 
diatomic molecules, known to exhibit this problem when studied within standard approaches, includ¬ 
ing the local spiir-density approximation (LSDA). By employing our recently proposed [E. Kraisler 
and L. Kronik, Phys. Rev. Lett. 110 , 126403 (2013)] ensemble-generalization we find that asymp¬ 
totic fractional dissociation is eliminated in all systems examined, even if the underlying exchange- 
correlation (xc) is still the LSDA. Furthermore, as a result of the ensemble generalization procedure, 
the Kohn-Sham potential develops a spatial step between the dissociated atoms, reflecting the emer¬ 
gence of the derivative discontinuity in the xc energy functional. This step, predicted in the past 
for the exact Kohn-Sham potential and observed in some of its more advanced approximate forms, 
is a desired feature that prevents any fractional charge transfer between the system’s fragments. 
It is usually believed that simple xc approximations such as the LSDA cannot develop this step. 
Our findings show, however, that ensemble generalization to fractional electron densities automati¬ 
cally introduces the desired step even to the most simple approximate xc functionals and correctly 
predicts asymptotic integer dissociation. 


I. INTRODUCTION 


Density-functional theory (DFT) [in is a widely-used 
theoretical framework for studying the electronic proper¬ 
ties of matter. It is usually employed by mapping the 
original N interacting-electron system into an equivalent 
Kohn-Sham (KS) system of N non-interacting particles 
subject to a common effective potential, such that the 
ground-state electron density n(r) is reproduced. While 
exact in principle, this mapping involves the exchange- 
correlation (xc) density functional, Exc{n{v)\, whose ex¬ 
act form is unknown. Therefore, it is always approximate 
in practice. 

Present-day density-functional approximations (DFAs) 
already make it widely applicable to a variety of many- 
electron systems in physics, chemistry, and materials sci¬ 
ence |10l - ll5l| . However, there remain numerous challenges 
that common DFAs fail to meet. A significant prob¬ 
lem, which has both formal and practical implications, 
is the so-called problem of fractional dissociation. It is 
most easily demonstrated when considering a neutral di¬ 
atomic molecule, AB. Upon dissociation, the molecule 
must break up into two neutral atoms, A and B, with an 
integer number of electrons, and Ng, on each [^ . 
This observation is known as principle of integer prefer¬ 
ence [l^. This fundamental principle is not reproduced 
by many DFAs. Instead, one finds that a system of two 
well-separated atoms often reaches its energy minimum 
when the number of electrons on each of the atoms is 
fractional: q electrons on atom A and Ng — q elec¬ 

trons on atom B, with q € (—1,1) (ibl - f^ . 

The significance of this failure lies beyond the accu¬ 
rate description of dissociation in diatomic systems. It 
indicates that common DFAs may fail to describe charge 
transfer in molecules and materials (see, e.g., [26l - l^ and 
references therein). Furthermore, theoretical analysis of 


the problem of fractional dissociation allowed an under¬ 
standing of fundamental properties of many-electron sys¬ 
tems [i3,|3a,[^: If N is allowed to be fractional, the total 
energy, E{N), of a many-electron system must possess a 
piecewise-liuear depeudence on N. As a result 
i? 3 ,c[u(r)] may exhibit a discontiuuity iu its derivative at 
integer N. 

It has been shown [l3,113, El] that the fractional disso¬ 
ciation problem occurs in approximate DFAs due to the 
deviation of their energy curves from piecewise linearity: 
standard (semi-)local calculations for atoms show a con¬ 
vex behavior of the energy Ea{N^-\-q), where a £ {A, B}. 
As a result, the total energy of a well-separated diatomic 
molecule, which can be expressed as 

EA...B{q) = EAiK + q) + EBiK-q), ( 1 ) 


may reach an unphysical minimum at q ^ 0 (see Fig. [T] 
below). This failure has been discovered as early as 1982 
by Perdew et al. [ 13 : for an infinitely stretched LiH 
molecule treated within the local density approximation, 
and has since been found in various molecules with differ- 
eut DFAs [h 


i EH . Perdew et al. tin have shown 
that this would not have happened had the atomic en¬ 
ergy curves, and as a result the molecular energy curve, 
been exactly piecewise linear. Then, in terms of the total 
energy, the molecular energy curve EA...B{q) would have 
possessed a non-analytical minimum at g = 0, obeying 
the principle of integer preference. In terms of the KS 
potential, a ’plateau’ in the vicinity of one of the atoms 
would have emerged, shifting the KS potential and the 
energy levels associated with that atom [H, m, E^IEl. 
Because all these desired features are absent in the stan¬ 
dard local density approximation, spurious transfer of 
fractional charge is not precluded. 

Recently, we have proposed an approximate ensemble 
generalization of the Hartree-exchange-correlation (Hxc) 
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functional [s^ . In contrast to the usual approach, in the 
generalized approach one does not insert the ensemble 
density, which integrates to a fractional N, into the stan¬ 
dard form of the approximate Hxc functional chosen. In¬ 
stead, the Hxc energy is now evaluated as a weighted sum 
of Hxc energies obtained using two auxiliary densities, 
which integrate to integer TV’s. These are derived from 
the same self-consistent KS potential by choosing the 
highest, partially-occupied KS energy level as either fully 
occupied or completely unoccupied. We have shown [531 - 
[55j | that this ensemble generalization restores, to a large 
extent, the piecewise-linearity property in the energy and 
simultaneously introduces an appropriate derivative dis¬ 
continuity into the xc potential in a natural manner, 
even when the underlying xc functional is as simple as 
the local spin-density approximation (LSDA). All this is 
achieved while neither introducing empiricism nor chang¬ 
ing the underlying functional form. This generalization 
would appear to be of direct relevance to the question 
of asymptotic fractional dissociation. It is therefore of 
much interest to examine whether and how it affects this 
important problem. 

Here we employ the ensemble-generalized LSDA, de¬ 
noted eLSDA, to well-separated diatomic molecules. We 
find that spurious asymptotic fractional dissociation is 
eliminated in all systems examined. Furthermore, the 
KS potential develops the desired spatial step, which ex¬ 
plains the absence of spurious fractional charge transfer 
also from the potential perspective. This shows that, 
contrary to conventional wisdom, the asymptotic frac¬ 
tional dissociation problem can be eliminated even with 
simple xc functionals, as long as an appropriate ensemble 
generalization is used. 


II. METHODOLOGY 

In the dissociation limit, a diatomic molecule can be 
described as two atoms with a possibly fractional elec¬ 
trical charge on each (Eq. ([T|)) [131 ■ Generally, a closed 
system with a fractional number of electrons, N = No+a, 
where iVo € N and a G ( 0 , 1 ) is not physical, in the sense 
that electrons do not fracture in chemical systems. A 
fractional electron number may, however, arise as a time 
average of an open system, which is free to exchange elec¬ 
trons with its surroundings. Such a system can no longer 
be described by a pure quantum-mechanical state. In¬ 
stead, it must be considered as a statistical mixture, or 
ensemble, of pure (integer electron) states [i3- Specifi¬ 
cally, in our case - Coulomb systems at zero temperature 
- this ensemble state is a linear combination of the pure 
ground states for Nq and A^o + 1 electrons, with the sta¬ 
tistical weights of (1 — a) and a, respectively [ 1 , 

[13, [H . 

For describing a quantum system with fractional IV 
in KS-DFT, recall that in KS-DFT one maps the origi¬ 
nal many-electron system to a fictitious non-interacting 
one, such that the overall density is retained. In par¬ 


ticular, the number of particles in the KS system must 
equal the number of electrons in the real system. There¬ 
fore, the KS density must also integrate to a fractional 
number of particles and the KS ground state must be 
an ensemble of Nq- and Nq + 1-states, i.e., Aj^s = 
(l-a)|<))(<)|+a|<Vi)(<Vil- Both |$M) and 
are Slater determinants constructed from one- 

electron KS orbitals, arising from the same KS 

potential. The difference between the two determinants 
is that is constructed from Nq orbitals, whereas 

is constructed from Nq + I orbitals. 

The above-discussed mapping is enforced by the KS 
potential, which is the common potential “felt” by all 
KS particles. Therefore, when changing the number of 
electrons, i.e., changing a, the KS potential itself can 
change because the nature of the mapping is changing. 
Significant variation of the KS potential with a have been 
indeed found for both the exact KS potential and various 
approximations to it (see [53 and references therein). As 
a result, the KS orbitals, v3-“^(r), and all the quantities 
derived from them (including (r) and (r) that are 
introduced below), are generally a-dependent. Here we 
emphasize this by including the superscript (a) explicitly. 

The usual approach to the KS-DFT treatment of sys¬ 
tems with fractional N is based on the construction of 
the electron density as 

No 

n(“)(r) = ^ |<p(“)(r)|2 + a|<pS+i(r)^ (2) 

i=l 

i.e., on occupying the hrst Nq levels fully but occupying 
the next one by the electron fraction a, so as to obtain 
No + a electrons overall. This density can be equivalently 
expressed as n(“^(r) = (1 — a)pQ°‘\r) + ap[‘^\r), where 
(r) is defined as the density due to the first Nq + p 
KS orbitals: Pp“^(r) = where p = 0 

or 1. In other words, the overall density is the weighted 
sum of the two auxiliary densities obtained from taking 
the partially occupied orbital as either completely empty 
or completely full (^ . These two auxiliary densities 
are not to be confused with the true densities of the Nq 
and Nq + 1 electron systems. In the usual approach, 
the overall density, n(“)(r), is then employed directly in 
the evaluation of the approximate Hxc energy functional, 
EhxcVA^ by inserting n*^“^(r) into the the functional form 
used for the description of systems with integer TV. 

Recently, we have shown [H, [13 that this usual ap¬ 
proach, i.e., using the integer Hxc form for systems with 
fractional TV, contributes significantly to the spurious de¬ 
viation from piecewise linearity in the total energy dis¬ 
cussed above. To address this problem, we proposed an 
approximate ensemble-based generalization of the Hxc 
functional to fractional TV. For completeness, we provide 
a short overview of this generalization. Its basic idea is 
that by considering the ensemble average of the many- 
electron Coulomb operator, W = 5 k* “ 
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in the KS system, the Hxc energy can be generalized as 
follows: 

= (1 - a)EHo.c[pi"^] + aEH.,[p[“\ (3) 

Here, the index e— signifies that the functional is 
ensemble-generalized, Ehxc is the pure-state Hxc func¬ 
tional mentioned above. The essential difference be¬ 
tween the ensemble-generalized functional Ef,- 
and the usual pure-state EHxc[n^°‘'^] is that we do not lin¬ 
early combine the pure state densities p^\r) hrst, and 
then inserting the result in the Hxc functional. Instead, 
following the ensemble approach, we first evaluate the 
Hxc for the pure-state densities separately, then linearly 
combine the two ensuing Ehxc[p^\tc)] energies. These 
two procedures are not the same mathematically, be¬ 
cause the Hxc functional is not linear with respect to 
the density. Only at integer N does the Hxc energy re¬ 
duce to the same form, namely that of the underlying 
pure-state Hxc functional. Because at fractional N the 
Hxc functional is explicitly dependent on /9p“^(r), i.e., it 
is explicitly orbital-dependent, treatment within the op¬ 
timized effective potential (OEP) formalism [3, IH, llfll 
is required. This is true even if the underlying xc func¬ 
tional is explicitly density-dependent, as is the case with 
the LSDA. 

The generalization in Eq. is applicable to any func¬ 
tional and makes the total energy explicitly linear in 
a. However, there may still remain an implicit non¬ 
linear dependence of the energy on a because of the 
a-dependent KS orbitals. In practice, the ensemble- 
generalization of Eq. ([3]) has already been shown to signifi¬ 
cantly improve, but not completely eliminate, deviations 
from piecewise linearity in the energy of simple atoms 
and molecules (s^j . 

Here we focus on infinitely separated diatomic 
molecules, which can be constructed from the atoms H, 
Li, C, and E, namely H...H, Li...H, C...H, F...H, Li...Li, 
Li...E, C...C, C...E and F...E, as well as their ions. All 
calculations are performed within the LSDA and eLSDA, 
while varying the number of electrons on each of the con¬ 
stituent atoms, using the DARSEC code [^. This code 
allows for spin-polarized all-electron DFT calculations 
for single atoms and diatomic molecules, using a real- 
space prolate-spheroidal grid. In all calculations the to¬ 
tal energy and the highest occupied (ho) eigenvalue have 
been obtained with an accuracy of 0.001 Ry. For all 
fractional densities, eLSDA calculations were performed 
within the optimized effective potential (OEP) formal¬ 
ism [3, m, 1^, usually in the Krieger-Li-Iafrate (KLI) 
approximation [^ . A special case that required em¬ 
ployment of the S-iteration method [o^, for solving 
the full OEP equation is discussed in the Supplementary 
Material [s^ . LSDA results are presented only if the ob¬ 
tained system is bound, which is not always the case for 
(fractionally) negatively charged ions. For eLSDA, where 
negative ions were never found to be bound, their energy 
is taken to be equal to that of a neutral system (plus 


a fraction of an electron at infinity, whose contribution 
is zero), [4l|, [M - as elaborated in the Supplementary 
Material [83] . 


III. RESULTS AND DISCUSSION 

We hrst consider the total energy of neutral molecules, 
Ea...b{q)-i as a function of g - the amount of negative 
charge transferred from B to A. Figure [1] presents the 
total energy of the Li and F atoms as a function of N 
and the total energy curve of the dissociated molecule 
Li...F as a function of q. The latter curve is obtained 
from the combination of the atomic curves, according to 
Eq. ([l|). For the atoms, a clear improvement in restoring 
piecewise linearity with eLSDA is observed, comparing 
to LSDA. For the dissociated molecule Li...F, the LSDA 
produces a convex energy curve, which is far from the 
expected piecewise-linear behavior, and with a spurious 
minimum at go = —0.4, in agreement with Ref. |l8| . In 
contrast, the eLSDA yields a result that is much closer 
to the piecewise-linear one, being a combination of two 
slightly concave curves. It possesses no spurious mini¬ 
mum, but rather a non-analytical minimum at 0, as re¬ 
quired. Similar spurious minima are obtained with LSDA 
for the neutral Li...H, F...H, and C...F molecules, with 
go = —0.05, 0.12, and —0.20, respectively, all removed 
by using eLSDA. These additional cases are shown in 
the Supplementary Material [83] . Importantly, close 
inspection of Fig. [T] (and its analogues in the Supple¬ 
mentary Material) reveals that the energy curves ob¬ 
tained with eLSDA are not perfectly linear but rather 
are slightly convex. This is a consequence of the remain¬ 
ing implicit a-dependence of the e-Hxc functional, owing 
to the above-discussed a-dependence of the KS orbitals, 
Also for ionized molecules of the type (A...B) + , 
where g = 0 corresponds to the state A“''...B and g = I 
to A...B+, with LSDA we find (C...H)+, (F...H)+, and 
(C...F)+ to possess spurious minima at go = 0.43, 0.67, 
and 0.27 respectively. With eLSDA all these minima are 
eliminated. For the homoatomic molecules H...H, Li...Li, 
C...C and F...F, as well as for C...H, (Li...H)+, (Li...F) + , 
no spurious minima occur even with the LSDA, but again 
the eLSDA energy curve is much closer to piecewise- 
linearity. 

Of particular interest are the E{q) curves for ho¬ 
moatomic ionized molecules (A...A)+ [13, For 

these systems we expect the energy to be the sum of the 
energies of the atom A and the ion A'*', independently 
of g. However, here the LSDA incorrectly prefers the 
state with g = 0.5, where the electron density is strongly 
delocalized, i.e., it is present on both nuclei. This is 
demonstrated, using (H...H)+ and (F...F)+ as examples, 
in Fig. [21 with more examples given in the Supplementary 
Material [83] . With LSDA, the energy for the g = 0.5 
state is too low by 0.20 and 0.29 Ry, respectively. The 
eLSDA cures this problem completely for (H...H)+, pro¬ 
ducing a flat energy curve, as required. For (F...F)+, the 
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FIG. 1: (Color online) Total energy of the Li and F atoms 
as a function of N (top) and energy of the dissociated Li...F 
molecule as a function of q (bottom), computed with LSD A 
(dash-dotted line) and eLSDA (solid line), compared to the 
expected piecewise-linear (dotted line) behavior (the latter 
is obtained by linear interpolation of the eLSDA energies at 
integer electron values). 


eLSDA results in a somewhat concave curve, which has 
its minima at q = 0 and 1, being too high by 0.1 Ry at 
q = 0.5. Thus, integer preference is again enforced. In 
fact it is ’’over-enforced”, in the sense that the integral 
states at q =0 and 1 are preferred over those with frac¬ 
tional q, instead of being iso-energetic with them. This 
is a consequence of the remaining deviation of the energy 
curve from a horizontal line, which is again traced to 
the implicit a-dependence of the e-Hxc functional. Still, 
spurious delocalization is completely eliminated. 

Fractional dissociation and its elimination can be fur¬ 
ther understood by examining dEA...B/dq. From Eq. (IT|) 
it follows that dEA...B/dq = ^J,A{q) — fJ-B{q), where /Tq := 
dEa/dN is the chemical potential of the a-th atom. From 
the property of piecewise-linearity of the exact E{N) we 
expect ^A “Ms to be a stair-step function of q. The de¬ 
pendence of Ha — Hb on q for Li...F is shown in Fig. [S] 
With LSDA we find that for low enough values of q, 
HA < Hb- This suggests that charge has to be moved 
from B to A, thereby increasing q. By doing so, we reach 
a chemical equilibrium for qo = —0.4, with ha = Hb [i3- 
For higher values oi q, ha > ^^B^ which suggests moving 
charge back from A to B, returning to the equilibrium 
at go- For eLSDA we find a completely different behav¬ 




FIG. 2: (Color online) Total energy of the (H...H)'*’ and 
(F...F)''’ molecules as a function of q computed with the LSDA 
(dash-dotted line) and eLSDA (solid line) functionals, com¬ 
pared to the expected piecewise-linear (dotted line) behavior. 


ior; ha{<i) — Hb{<i) is negative for negative g’s, suggesting 
a transfer of charge from B to A, i.e., increasing g, until 
reaching 0, where a discontinuous change is observed and 
Ha{<i) — ^^B{^) becomes positive, suggesting decreasing g 
back. This illustrates the fact that E{q) exhibits a non- 
analytic minimum: while the system reaches its energy 
minimum at g = 0, the two atoms are, strictly speaking, 
never in chemical equilibrium, because ha ^ Hb- An in¬ 
finitesimal increase of g above g = 0 results in an abrupt 
jump, which indicates that g has to be decreased, and 
vice versa. Note that the deviation of the eLSDA from 
the expected stepwise form in Fig. [3] originates from the 
residual convexity of the E{q) in Fig. [TJ 

The energy considerations presented above are also re¬ 
flected in the KS potential. For an infinitely separated 
diatomic molecule, the molecular potential coincides with 
the atomic potential in the vicinity of each of the nuclei. 
Therefore, here it is obtained by combination of the two 
KS atomic potentials obtained from two separate atomic 
calculations. 

As a result of the ensemble generalization to the energy 
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FIG. 3: (Color online) Dependence of fJ^A{q) — Ms(9) on 9 for 
the Li...F molecule, computed with LSDA (dash-dotted line) 
and eLSDA (solid line), compared to the expected discontin¬ 
uous stair-step (dotted line) behavior. 


(Eq. dS])), the KS potential is generalized as well. Tak¬ 
ing the functional derivative of Ee-Hxc[n\ with respect 
to n, while recalling that the ensemble-generalized func¬ 
tional is both a- and orbital-dependent one reveals 
a peculiar property of the KS potential 53|: it does not 
necessarily tend asymptotically to 0, but rather to a fi¬ 
nite spatial constant, Voo- We stress that this constant is 
well-dehned and can be analytically expressed in terms 
of the KS orbitals and Hxc potential [5J, [dj . It must be 
taken into account in order for the ho KS eigenvalue to 
equal dE/dN, i.e., to obey Janak’s theorem 1^. Indeed, 
including this constant signihcantly improves |53l. Issjl the 
prediction of the ionization potentia l (IP) from the ho 
KS eigenvalue via the IP theorem [13, ItoI - I^ . In the cur¬ 
rent context, we emphasize that Uoo is system-dependent 
and may change discontinuously with N across an integer 
point. 

To illustrate how the ensemble generalization prevents 
spurious charge transfer from the potential perspective, 
we use the (C...F)++ molecule. We prefer it to the Li...F 
used above to confine ourselves to a system with strictly 
bound states (see the Supplementary Material [83] for 
a complete discussion). 

The energy curve of the system is qualita¬ 

tively similar to the one shown at Fig. [U with a spurious 
minimum at qo = —0.273 for LSDA and a non-analytic 
minimum at 0 for eLSDA. In the following, we discuss 
the quantity w(r), which equals the J,-KS potential, mi¬ 
nus the external and the Hartree potentials. The latter 
do not possess any discontinuity and therefore are sub¬ 
tracted for clarity. We stress that Fig. 0] that depicts 
w{r) is not at all a schematic drawing. It is obtained by 
combining two KS potentials obtained from two separate 
atomic calculations, for and F^+^, while preserving 
their asymptotic values [s^j . 



FIG. 4: (Color online) The potential w{r) for the system 
C^“'^...F^"'"’ plotted against the interatomic axis, z: (a) with 
LSDA, at q — qo = —0.273; (b) with eLSDA, at q = —0.001; 
(c) with eLSDA, at q = 4-0.001. The ho levels of both ions 
are marked with horizontal lines. 


For LSDA (Fig. |4][a)), w{r) essentially differs from 0 
only near the two nuclei and decays to 0 elsewhere. For 
q < qo, the ho level of the C ion is lower than the ho 
level of the F ion. The levels approach each other gradu¬ 
ally as q approaches qo, indicating chemical equilibrium. 
Further increase of q makes the levels gradually sepa¬ 
rate. Therefore, chemical equilibrium is achieved by a 
spurious transfer of a fraction of an electron from the 
C ion to the F ion, contrary to the principle of integer 
preference. For eLSDA (Figs. HKb) and (c)), the situa¬ 
tion is radically different. Since the ensemble-generalized 
KS potential tends asymptotically to a non-zero system- 
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dependent constant, Voo-, the potential experiences a step 
between the two well-separated ions. Moreover, since Voo 
is discontinuous when N passes through an integer, the 
value of the step will abruptly change around an integer 
q. For q < 0, we observe a positive step, s_, so the ho of 
the C ion is always lower than that of the F ion, suggest¬ 
ing that q should be increased. However, for g > 0, the 
step s+ is negative and the ho of the C ion is higher than 
that of the F ion, suggesting that q should be decreased. 
This mechanism of negative reaction to charge transfer 
assures that the well-separated subsystems will always be 
occupied by an integer number of electrons, as expected. 
The height of the steps can be expressed analytically in 
KS quantities of the subsystems A and B with inte gra l 
N: s- = v% — v^ + Ab; s+ = v% — — A^, where 

Va — E!Hxc[na] ~ ^ | ] 

- J d^r\ip^hlai^)\^^Hxc[na]{r) (4) 

is the ensemble-KS-potential asymptote for Na (A"a)~ 
(i.e., a 1“) [13, and 

Aa = EHxc[naM‘p\l]a\‘^]-‘^EHxc[na]+EHxc[na-\‘Phla\'^] 

+ J d^r UB,,cK](r) (5) 

is the ensemble derivative discontinuity of the a-th sys¬ 
tem (a S {A,B}) that occurs as N crosses the integer 
value Note that for the exact functional = 0, i.e., 
the asymptotic potentials tends to zero, and therefore 
s- = Ab and s+ = —A^- 

Importantly, we note that here we only addressed in¬ 
finitely separated molecules, by treating their constituent 
atoms independently. Naive application of the above- 
formulated ensemble generalization to finite bond-length 
molecules is not useful as the overall number of electrons 
in the system is still integer, instigating no ensemble cor¬ 
rection. Further generalization of the approach, possibly 
by employing partition-DFT to address sub-systems (see, 
e.g., and references therein), so as to overcome this 
difficulty will be discussed elsewhere. 


IV. CONCLUSIONS 

In conclusion, from an energy perspective the problem 
of asymptotic fractional dissociation in DFT emerges ow¬ 
ing to significant convexity in the energy versus fractional 
electron number curves. Here, we have shown that this 
convexity is replaced by much smaller concavity by em¬ 
ploying the ensemble generalization (s^ . thereby restor¬ 
ing integer dissociation, even when the underlying xc 
functional is the simple LSDA. 

From a potential perspective, spurious fractional 
charge transfer between well-separated molecular frag¬ 
ments occurs in standard approximations because they 
lack a ’plateau’ in the KS potential. This ‘plateau’ is nat¬ 
urally introduced by the ensemble generalization. There¬ 
fore, the ensemble generalization enforces the physical 
integer dissociation limit in both the energy and the po¬ 
tential pictures, without introducing empiricism or al¬ 
tering the underlying functional form. Finally, we note 
that integer dissociation is recovered by this procedure 
despite the fact that both LSDA and eLSDA suffer from 
one-electron self-interaction. We therefore conclude that 
the connection between fractional dissociation and one- 
electron self-interaction [Il,[zi,[z3 is not a straightfor¬ 
ward one. In our opinion, self-interaction correction 
schemes can indeed remedy fractional dissociation, but 
this is because they produce a concave energy curve. We 
expect that the ensemble-generalization perspective can 
bring significant improvement in the description of charge 
transfer processes, even with currently existing xc func¬ 
tionals. 
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I. ENERGY CURVES E{N) FOR THE ATOMS H, Li, C, AND F 

Figure [T] presents the dependence of the energy, E{N), on the number of electrons, N, obtained for various atoms, 
using the LSDA and eLSDA. The results are compared to the piecewise-linear behavior obtained by linear interpolation 
between eLSDA values at integer N. Data from these curves were used to obtain the energy curves of infinitely 
separated molecules, as shown in Sec. m below. 


0 

- 0.2 
-0.4 
£ - 0.6 
- 0.8 
-1 

- 1.2 



N 


N 


-72.5 


-73.5 

Pi 




-74.5 


-75.5 



FIG. 1: Energy as a function of the number of electrons, for the atoms H, Li, C, and F, computed with LSDA (dash-dotted 
line) and eLSDA (solid line), compared to the expected piecewise-linear (dotted line) behavior (the latter is obtained by linear 
interpolation of the eLSDA energies at integer electron values). 
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From Fig. [T] one can clearly see that the eLSDA is much closer to piecewise-linearity in all cases. As mentioned 
in the main text, LSDA results are presented only if the obtained system is bound, which is not always the case for 
(fractionally) negatively charged ions. For eLSDA, where negative ions were never obtained as bound, their energy 
was taken to be equal to that of a neutral system (plus a fraction of an electron at infinity, whose contribution is 
zero) An alternative approach to negative ion calculations is discussed in Sec. [TTIl below. Note that in the 

context of the fractional dissociation problem both approaches lead to the same conclusions. 


II. ENERGY CURVES E{q) FOR INFINITELY SEPARATED MOLECULES 

Using the data from atomic calculations, presented in Fig. [U we can obtain energy curves Ea...b((i) for infinitely 
separated molecules, using Eq. (1) from the main text, q is the amount of negative charge transferred from atom B 
to atom A. Figure [2] presents the energy curves for neutral, heteroatomic molecules. Spurious mimina obtained with 
the LSDA are marked by rhombi. We readily see that these minima are absent from eLSDA calculations in all cases 
examined. Figure [3] presents the energy curves for neutral, homoatomic molecules. For these systems, no spurious 
minima are expected, due to symmetry considerations. However, we see that eLSDA is much closer to piecewise- 
linearity. Figured presents the energy curves for singly ionized, heteroatomic molecules and Fig. [5] for singly ionized, 
homoatomic molecules. Spurious mimina obtained with the LSDA are marked by rhombi. We again readily see that 
these minima are absent from eLSDA calculations in all cases examined. Finally, Fig. [S] shown the energy curve for 
the (C...F)++ molecule, which has been chosen in the main text to study the KS potential behavior. 




FIG. 2: Energy as a function of q for the neutral, heteroatomic molecules Li...H, C...H, F...H, and C...F, computed with 
LSDA (dash-dotted line) and eLSDA (solid line), compared to the expected piecewise-linear (dotted line) behavior (the latter 
is obtained by linear interpolation of the eLSDA energies at integer electron values). 
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FIG. 3: Energy as a function of q for the neutral, homoatomic molecules Li...Li, C...C, and F...F, computed with LSDA 

(dash-dotted line) and eLSDA (solid line), compared to the expected piecewise-linear (dotted line) behavior (the latter is 
obtained by linear interpolation of the eLSDA energies at integer electron values). 
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FIG. 4: Energy as a function of q for the singly ionized, heteroatomic molecules (Li...H)'^, (Li...F)''", and 

(C...F)'^, computed with LSDA (dash-dotted line) and eLSDA (solid line), compared to the expected piecewise-linear (dotted 
line) behavior (the latter is obtained by linear interpolation of the eLSDA energies at integer electron values). 
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FIG. 5: Energy as a function of q for the singly ionized, homoatomic molecules (Li...Li)'^ and (C...C)’^, computed with LSDA 
(dash-dotted line) and eLSDA (solid line), compared to the expected piecewise-linear (dotted line) behavior (the latter is 
obtained by linear interpolation of the eLSDA energies at integer electron values). 



FIG. 6: Energy as a function of q for the doubly ionized molecule (C...F)’'"^, computed with LSDA (dash-dotted line) and eLSDA 
(solid line), compared to the expected piecewise-linear (dotted line) behavior (the latter is obtained by linear interpolation of 
the eLSDA energies at integer electron values). 


III. APPROACHES TO NEGATIVELY CHARGED IONS 

In order to describe an infinitely separated diatomic molecule, while transferring a fraction of an electron from one 
of its atom to the other one, we sometimes need to perform atomic calculations for ions with a negative net charge, 
i.e., when the (possibly fractional) number of electrons, N, is larger than the atomic number, Z. It is well-known that 
the LSDA and eLSDA may yield unbound solutions for such systems Q. 

Principally, there are two approaches to address the aforementioned case ot N > Z (i) To state that since the 

negative ion is not bound, its energy, as well as the energy of an ion with any fractional N larger than Z, equals the 
energy of the neutral atom plus the energy of (a fraction of) an electron at infinity. Since the contribution of the latter 
is zero, for iV > Z the energy E{N) = E{Z), i.e., it is a horizontal line, independent of N; (ii) To obtain the energy 
of the negatively charged system from a DFT calculation (see, e.g., and references therein). Here, however, we 
have to note that in cases where the system is not bound, calculations are normally not converged in terms of the 
basis set, or in our case of real-space calculations - in terms of fj-max ~ the maximal value of the prolate-spheroidal 
coordinate /i, which plays here roughly the same role as the radial coordinate r in spherical coordinates. Therefore, 
usually the grid of the neutral system is used for negative ions, as well, although it is formally not converged. 






















6 


In the current work, we chose approach (i). In this section we show that employment of approach (ii) to total- 
energy calculations of infinitely-separated molecules with LSDA and eLSDA, although yeilding results that are different 
quantitatively, leads to the same qualitative conclusions as approach (i): removal of the spurious energy minimum 
at fractional q by eLSDA and strongly reduced deviation from piecewise linearity. For illustration we consider the F 
atom and the molecule Li...F. 

In Fig. Owe compare the E{N) curves of the F atom, obtained with LSDA and eLSDA, for both approaches (i) and 
(ii). For LSDA, we find that for N > 9.8 an unbound solution is obtained with approach (ii). For eLSDA, approach 
(ii) yields a convex energy curve, which approaches the same energy value as the LSDA curve at A^ = 10. For the 
other atoms, namely H, Li, and C, a similar picture is obtained. 



FIG. 7: Energy as a function of N for the F atom, obtained using: approach (i) - with LSDA (blue, closed squares, solid line) 
and eLSDA (red, closed rhombuses, solid line), compared to the expected piecewise-linear (black, closed circles) behavior (the 
latter is obtained by linear interpolation of the eLSDA energies at integer electron values); approach (ii) - with LSDA (blue, 
open squares, dashed line) and eLSDA (red, open rhombs, dashed line), compared to the expected piecewise-linear (black, open 
circles, dashed line) behavior (the latter is obtained by linear interpolation of the eLSDA energies at integer electron values). 

In Fig.|S]we present the E{q) curves for the Li...F molecules, obtained with both approaches. We clearly see that the 
spurious LSDA energy minimum at qo = —0.4 is eliminated by the eLSDA calculations performed in both approaches 
(i) and (ii). Instead, they produce a non-analytical minimum at q = 0, as required. 



FIG. 8: Energy as a function of q for the Li...F molecule, obtained using approaches (i) and (ii). All symbols are as in Fig. [T] 
The piecewise-linear behavior in approach (ii) is plotted only for q < 0, for clarity 

We note that employment of approach (ii) in the calculations included several significant numerical difficulties. 
First, for the H and Li negative ions we found that if using DARSEC’s standard initial guess for the density, for a 
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certain N the total energy jumps abruptly upwards and the energy levels drop downwards, yielding an artificial 
solution that is bound on one hand, but is higher in energy than that for the neutral atom. For Li, this happened at 
Ni = 3.346. We relate this phenomenon to the fact that the numerical grid used in the calculation is not converged. 
To overcome the appearance of this solution, we had to restart our calculations with an initial guess taken from a 
calculation with N < Ni. 

Second, for the C and F negative ions, the self-consistent cycle could not converge when using the KLI approxima¬ 
tion, due to the proximity of the p-levels in these atoms and frquent ch ang es in their order during the self-consistent 
cycle. To overcome this difficulty, we had to use the S-iteration method (7|, l8| . Since the latter method is numerically 
sensitive, it is usually employed only after performing an initial step, e.g., of KLI. In our case, we used the simple 
Slater potential [1,0 for the first iterations. 

To summarize, a principal deficiency of the LSDA and eLSDA, i.e., not binding negative ions, raises numerous 
technical difficulties that obscures the discussion of the fractional dissociation problem. Ways to overcome them have 
been discussed above. Note that while relevant for neutral well-separated molecules, calculations of negative ions are 
not required for positively charged molecules. The latter have been considered in our work especially for illuminating 
the fractional dissociation problem in systems where binding problems are absent. As shown above, in positively 
charged molecule the fractional dissociation problem is removed by the ensemble treatment just as well as in neutral 
systems. 
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